home *** CD-ROM | disk | FTP | other *** search
/ Atari Mega Archive 1 / Atari Mega Archive - Volume 1.iso / apps / astronmy / strchart.arc / patch.arc / planet.dif < prev    next >
Text File  |  1989-03-23  |  18KB  |  598 lines
  1. *** ../starchart/planet.c    Mon Sep 19 09:11:43 1988
  2. --- planet.c    Mon Sep 26 08:00:31 1988
  3. ***************
  4. *** 45,57 ****
  5.   #include <stdio.h>
  6.   #include <math.h>
  7.   
  8.   #ifndef SYSV
  9.   #include <sys/time.h>    /* for getting current GMT (generic Unix) */
  10.   #else
  11.   #include <time.h>    /* for getting current GMT (sysV version) */
  12.   #endif
  13.   
  14. ! #define CURYEAR    1987    /* default year -- needs to be maintained */
  15.   
  16.   #ifndef PLANETFILE
  17.   #define PLANETFILE "./planet.star"
  18. --- 45,62 ----
  19.   #include <stdio.h>
  20.   #include <math.h>
  21.   
  22. + #ifdef ST_MWC /* Mark Williams C for the Atari ST */
  23. + #define SYSV /* MWC is almost SYSV compatible */
  24. + #include <osbind.h> /* Needed for Cconws() and Crawcin() */
  25. + #endif /* ST_MWC */
  26. + 
  27.   #ifndef SYSV
  28.   #include <sys/time.h>    /* for getting current GMT (generic Unix) */
  29.   #else
  30.   #include <time.h>    /* for getting current GMT (sysV version) */
  31.   #endif
  32.   
  33. ! #define CURYEAR    1989    /* default year -- needs to be maintained */
  34.   
  35.   #ifndef PLANETFILE
  36.   #define PLANETFILE "./planet.star"
  37. ***************
  38. *** 68,73 ****
  39. --- 73,84 ----
  40.   #define MAGURA 590
  41.   #define MAGNEP 800
  42.   
  43. + #ifdef ST_MWC
  44. + #define daylight tt.tm_isdst
  45. + extern char *getenv(); /* used to see if we're running from the desktop */
  46. + struct tm_t *tt;
  47. + #endif /* ST_MWC */
  48. + 
  49.   double pie, rad;
  50.   double htod(), atof(), kepler(), truean();
  51.   double longi(), lati(), poly(), aint(), range();
  52. ***************
  53. *** 119,125 ****
  54. --- 130,138 ----
  55.       gettimeofday(&tv, &tz);
  56.       return(GMT1970 + tv.tv_sec/SECSPERDAY);
  57.   #else
  58. + #ifndef ST_MWC
  59.       long time();
  60. + #endif /* ST_MWC */
  61.       return(GMT1970 + (double)time((long *)0)/SECSPERDAY);
  62.   #endif
  63.       }
  64. ***************
  65. *** 201,207 ****
  66.       aa1 = (int) (year/100);
  67.       bb1 = 2 - aa1 + (int)(aa1/4);
  68.       jd = aint(365.25*year) + aint(30.6001*(month + 1));
  69. !     jd = jd  + day + 1720994.5;
  70.       if((year + month/100) > 1582.10) jd = jd + bb1;
  71.       return(jd);
  72.       }
  73. --- 214,220 ----
  74.       aa1 = (int) (year/100);
  75.       bb1 = 2 - aa1 + (int)(aa1/4);
  76.       jd = aint(365.25*year) + aint(30.6001*(month + 1));
  77. !     jd += day + 1720994.5;
  78.       if((year + month/100) > 1582.10) jd = jd + bb1;
  79.       return(jd);
  80.       }
  81. ***************
  82. *** 218,224 ****
  83. --- 231,243 ----
  84.       double N,D,epli,thapp,omeg;
  85.       double nu2,P,Q,S,V,W,ze,l1pert,epert,w1pert,apert;
  86.       double psi,H,G,eta,th;
  87. + #ifdef ST_MWC
  88. +     time_t now;
  89.   
  90. +     now = time((long) 0);
  91. +     tt = localtime(&now);
  92. + #endif /* ST_MWC */
  93. + 
  94.   #define WRITEMODE "w"
  95.   #define OPENFAIL 0
  96.       progname = argv[0];
  97. ***************
  98. *** 321,327 ****
  99.           RA = atan2(N,D)/rad;
  100.           DEC = asin(sin(epli*rad)*sin(thapp*rad))/rad;
  101.           speak(RA,DEC,Sr, MAGSOL, "PS", "Sol");
  102. ! /* tansformation of coordinates on Mercury and output */
  103.           trans(r,b,ll,Stheta,Sr,epli, MAGMER, "PM", "Mercury");
  104.           
  105.   /* Now start on Venus */
  106. --- 340,346 ----
  107.           RA = atan2(N,D)/rad;
  108.           DEC = asin(sin(epli*rad)*sin(thapp*rad))/rad;
  109.           speak(RA,DEC,Sr, MAGSOL, "PS", "Sol");
  110. ! /* transformation of coordinates on Mercury and output */
  111.           trans(r,b,ll,Stheta,Sr,epli, MAGMER, "PM", "Mercury");
  112.           
  113.   /* Now start on Venus */
  114. ***************
  115. *** 373,380 ****
  116.       radpert = 0.000022501 * cos((2*M-2*M2-58.208)*rad)
  117.            +0.000019045 * cos((3*M-3*M2+92.577)*rad)
  118.            +0.000006887 * cos((M5-M2-118.090)*rad)
  119. !          +0.000005172 * cos((M-M2-29.110)*rad)
  120. !          +0.000003620 * cos((5*M-4*M2-104.208)*rad)
  121.            +0.000003283 * cos((4*M-4*M2+63.513)*rad)
  122.            +0.000003074 * cos((2*M5-2*M2-55.167)*rad);
  123.       r = r + radpert;
  124. --- 392,399 ----
  125.       radpert = 0.000022501 * cos((2*M-2*M2-58.208)*rad)
  126.            +0.000019045 * cos((3*M-3*M2+92.577)*rad)
  127.            +0.000006887 * cos((M5-M2-118.090)*rad)
  128. !          +0.000005172 * cos((M-M2-29.110)*rad);
  129. !     radpert+= 0.000003620 * cos((5*M-4*M2-104.208)*rad)
  130.            +0.000003283 * cos((4*M-4*M2+63.513)*rad)
  131.            +0.000003074 * cos((2*M5-2*M2-55.167)*rad);
  132.       r = r + radpert;
  133. ***************
  134. *** 425,432 ****
  135.       lonpert = 0.00705*cos((M5-M4-48.958)*rad)
  136.            +0.00607*cos((2*M5-M4-188.350)*rad)
  137.            +0.00445*cos((2*M5-2*M4-191.897)*rad)
  138. !          +0.00388*cos((M-2*M4+20.495)*rad)
  139. !          +0.00238*cos((M-M4+35.097)*rad)
  140.            +0.00204*cos((2*M-3*M4+158.638)*rad)
  141.            +0.00177*cos((3*M4-M2-57.602)*rad)
  142.            +0.00136*cos((2*M-4*M4+154.093)*rad)
  143. --- 444,451 ----
  144.       lonpert = 0.00705*cos((M5-M4-48.958)*rad)
  145.            +0.00607*cos((2*M5-M4-188.350)*rad)
  146.            +0.00445*cos((2*M5-2*M4-191.897)*rad)
  147. !          +0.00388*cos((M-2*M4+20.495)*rad);
  148. !     lonpert+= 0.00238*cos((M-M4+35.097)*rad)
  149.            +0.00204*cos((2*M-3*M4+158.638)*rad)
  150.            +0.00177*cos((3*M4-M2-57.602)*rad)
  151.            +0.00136*cos((2*M-4*M4+154.093)*rad)
  152. ***************
  153. *** 445,452 ****
  154.       radpert+= 0.000007914*cos((3*M5-2*M4-139.737)*rad)
  155.            +0.000007004*cos((2*M5-3*M4-102.888)*rad)
  156.            +0.000006620*cos((M-2*M4+113.202)*rad)
  157. !          +0.000004930*cos((3*M5-3*M4-76.243)*rad)
  158. !          +0.000004693*cos((3*M-5*M4+190.603)*rad)
  159.            +0.000004571*cos((2*M-4*M4+244.702)*rad)
  160.            +0.000004409*cos((3*M5-M4-115.828)*rad);
  161.       r = r + radpert;
  162. --- 464,471 ----
  163.       radpert+= 0.000007914*cos((3*M5-2*M4-139.737)*rad)
  164.            +0.000007004*cos((2*M5-3*M4-102.888)*rad)
  165.            +0.000006620*cos((M-2*M4+113.202)*rad)
  166. !          +0.000004930*cos((3*M5-3*M4-76.243)*rad);
  167. !     radpert+= 0.000004693*cos((3*M-5*M4+190.603)*rad)
  168.            +0.000004571*cos((2*M-4*M4+244.702)*rad)
  169.            +0.000004409*cos((3*M5-M4-115.828)*rad);
  170.       r = r + radpert;
  171. ***************
  172. *** 507,513 ****
  173.       -(0.003083 + 0.000275*nu2 - 0.000489*nu2*nu2)*sin(2*V)
  174.       +0.002472*sin(W)
  175.       +0.013619*sin(ze)
  176. !     +0.018472*sin(2*ze)
  177.       +0.006717*sin(3*ze)
  178.       +0.002775*sin(4*ze)
  179.       +(0.007275 - 0.001253*nu2)*sin(ze)*sin(Q)
  180. --- 526,533 ----
  181.       -(0.003083 + 0.000275*nu2 - 0.000489*nu2*nu2)*sin(2*V)
  182.       +0.002472*sin(W)
  183.       +0.013619*sin(ze)
  184. !     +0.018472*sin(2*ze);
  185. ! l1pert = l1pert
  186.       +0.006717*sin(3*ze)
  187.       +0.002775*sin(4*ze)
  188.       +(0.007275 - 0.001253*nu2)*sin(ze)*sin(Q)
  189. ***************
  190. *** 519,525 ****
  191.       -(0.035681 + 0.001208*nu2)*sin(ze)*cos(Q)
  192.       -0.004261*sin(2*ze)*cos(Q)
  193.       +0.002178*cos(Q)
  194. !     +(-0.006333 + 0.001161*nu2)*cos(ze)*cos(Q)
  195.       -0.006675*cos(2*ze)*cos(Q)
  196.       -0.002664*cos(3*ze)*cos(Q)
  197.       -0.002572*sin(ze)*sin(2*Q)
  198. --- 539,546 ----
  199.       -(0.035681 + 0.001208*nu2)*sin(ze)*cos(Q)
  200.       -0.004261*sin(2*ze)*cos(Q)
  201.       +0.002178*cos(Q)
  202. !     +(-0.006333 + 0.001161*nu2)*cos(ze)*cos(Q);
  203. ! l1pert = l1pert
  204.       -0.006675*cos(2*ze)*cos(Q)
  205.       -0.002664*cos(3*ze)*cos(Q)
  206.       -0.002572*sin(ze)*sin(2*Q)
  207. ***************
  208. *** 532,538 ****
  209.       -.0006764*sin(ze)*sin(Q)
  210.       -.0001110*sin(2*ze)*sin(Q)
  211.       -.0000224*sin(3*ze)*sin(Q)
  212. !     -.0000204*sin(Q)
  213.       +(.0001284 + .0000116*nu2)*cos(ze)*sin(Q)
  214.       +.0000188*cos(2*ze)*sin(Q)
  215.       +(.0001460 + .0000130*nu2)*sin(ze)*cos(Q)
  216. --- 553,560 ----
  217.       -.0006764*sin(ze)*sin(Q)
  218.       -.0001110*sin(2*ze)*sin(Q)
  219.       -.0000224*sin(3*ze)*sin(Q)
  220. !     -.0000204*sin(Q);
  221. ! epert = epert
  222.       +(.0001284 + .0000116*nu2)*cos(ze)*sin(Q)
  223.       +.0000188*cos(2*ze)*sin(Q)
  224.       +(.0001460 + .0000130*nu2)*sin(ze)*cos(Q)
  225. ***************
  226. *** 544,550 ****
  227.       +.0000508*cos(3*ze)*cos(Q)
  228.       +.0000230*cos(4*ze)*cos(Q)
  229.       +.0000108*cos(5*ze)*cos(Q)
  230. !     -(.0000956 + .0000073*nu2)*sin(ze)*sin(2*Q)
  231.       +.0000448*sin(2*ze)*sin(2*Q)
  232.       +.0000137*sin(3*ze)*sin(2*Q)
  233.       +(-.0000997 + .0000108*nu2)*cos(ze)*sin(2*Q)
  234. --- 566,573 ----
  235.       +.0000508*cos(3*ze)*cos(Q)
  236.       +.0000230*cos(4*ze)*cos(Q)
  237.       +.0000108*cos(5*ze)*cos(Q)
  238. !     -(.0000956 + .0000073*nu2)*sin(ze)*sin(2*Q);
  239. ! epert = epert
  240.       +.0000448*sin(2*ze)*sin(2*Q)
  241.       +.0000137*sin(3*ze)*sin(2*Q)
  242.       +(-.0000997 + .0000108*nu2)*cos(ze)*sin(2*Q)
  243. ***************
  244. *** 554,560 ****
  245.       +(-.0000956 +.0000099*nu2)*sin(ze)*cos(2*Q)
  246.       +.0000490*sin(2*ze)*cos(2*Q)
  247.       +.0000158*sin(3*ze)*cos(2*Q)
  248. !     +.0000179*cos(2*Q)
  249.       +(.0001024 + .0000075*nu2)*cos(ze)*cos(2*Q)
  250.       -.0000437*cos(2*ze)*cos(2*Q)
  251.       -.0000132*cos(3*ze)*cos(2*Q);
  252. --- 577,584 ----
  253.       +(-.0000956 +.0000099*nu2)*sin(ze)*cos(2*Q)
  254.       +.0000490*sin(2*ze)*cos(2*Q)
  255.       +.0000158*sin(3*ze)*cos(2*Q)
  256. !     +.0000179*cos(2*Q);
  257. ! epert = epert
  258.       +(.0001024 + .0000075*nu2)*cos(ze)*cos(2*Q)
  259.       -.0000437*cos(2*ze)*cos(2*Q)
  260.       -.0000132*cos(3*ze)*cos(2*Q);
  261. ***************
  262. *** 562,568 ****
  263.   w1pert = (0.007192 - 0.003147*nu2)*sin(V)
  264.       +(-0.020428 - 0.000675*nu2 + 0.000197*nu2*nu2)*cos(V)
  265.       +(0.007269 + 0.000672*nu2)*sin(ze)*sin(Q)
  266. !     -0